Improved position measurement of nanoelectromechanical systems using cross 
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We consider position measurements using the cross-correlated output of two tunnel junction posi- 
tion detectors. Using a fully quantum treatment, we calculate the equation of motion for the density 
matrix of the coupled detector-detector-mechanical oscillator system. After discussing the presence 
of a bound on the peak-to-background ratio in a position measurement using a single detector, we 
show how one can use detector cross correlations to overcome this bound. We analyze two different 
possible experimental realizations of the cross correlation measurement and show that in both cases 
the maximum cross-correlated output is obtained when using twin detectors and applying equal bias 
to each tunnel junction. Furthermore, we show how the double-detector setup can be exploited to 
drastically reduce the added displacement noise of the oscillator. 

PACS numbers: 85.85.+j,73.23.-b,72.70.-|-m 
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I. INTRODUCTION 

It is expected that, in the near future, position measure- 
ments of nanomechanical systems will reach the quantum 
Umit. Experimental progress in this direction is very fast 
and displacement sensitivities near the standard quan- 
tum limit have already been demonstrated^i^i^. 
In the current generation of experiments, the coupling be- 
tween the resonator and the mesoscopic detector is typi- 
cally very weak. The position measurement can therefore 
not be seen as a strong projective measurement. It is bet- 
ter described within the framework of weak measurement 
theory that was recently developed in the context of solid- 
state quantum computing^i^i^ii. This theory describes a 
continuous measurement process where the information 
about the measured object can be extracted, for instance, 
from the spectral density of the detector (and not simply 
from its average output). An important result in this the- 
ory is the Korotkov-Averin bound, which puts an upper 
limit of 4 to the ratio of the contribution of the measured 
state to the detector's spectral density, and the intrinsic 
background detector noise, for any linear detector mea- 
suring a two-level system. 

Since a quantum position measurement by a mesoscopic 
detector can be described within the same theoretical 
framework as a qubit measurement, one might ask if 
such a bound also exists in the case of a position mea- 
surement. In this article, we first show that, for fixed 
system parameters, the peak-to-background ratio in the 
spectral density of a position detector weakly coupled to 
an oscillator is also bounded from above. This result is 
obtained by considering the example of a single tunnel- 
junction detector, a simple detector that has been thor- 
oughly studied theoreticall y^'^i^°'^^i^^'^'^i^^ and realized 
exp er iment allyi^ii^ . 

Besides showing that the peak-to-background ratio is 
bounded in the typical single-detector position measure- 
ment, we also propose, in this article, two simple ex- 
perimental configurations (Fig. [1]) where, by using the 
cross correlations between two detectors, the bound on 



the peak-to-background ratio can be overcome. As the 
oscillator-independent parts of the output signal of the 
two detectors are uncorrelated, the background noise in 
these configurations is zero and therefore the peak-to- 
background ratio diverges. In the context of qubit read- 
out, this idea has already been proposed in an insightful 
work by Jordan and Biittiker— and was shown experi- 
mentally to improve readout fidelitj^. Experimentally, 
position measurements should hence also profit from us- 
ing cross-correlated detector outputs. We analyze in de- 
tail the two configurations presented in Fig. [T] and obtain 
analytical results for the optimal cross-correlated signal 
as a function of different detector parameters. 
Previous studies^^ of the position measurement problem 
focused on finding the conditions for quantum-limited de- 
tector sensitivity, under which one minimizes the total 
detector contribution to the output displacement noise. 
We show that the double-detector setup proposed here 
can in fact be used to almost totally get rid off the 
added displacement noise of the oscillator due to detec- 
tor back-action. This is a remarkable result that nicely 
com plem ents the general single-detector analysis made in 
Ref.lli. 

The article is organized as follows: in Section|TTl we intro- 
duce the formalism used in the rest of the paper, viz., a 
master equation for the m-resolved density matrix, where 
m is the number of charges that have passed through the 
detector. This equation of motion allows us to find ex- 
pressions for the combined moments of charge (detector) 
and oscillator quantities. In Section IIIIl the formalism 
is applied to the case of one position detector coupled to 
the oscillator. ''S, We analyze the peak-to-background ra- 
tio and show that this quantity is always bounded from 
above in the single-detector case. This bound cannot be 
made arbitrarily large simply by increasing the detector 
sensitivity. Section IIVI generalizes this treatment to a 
configuration with two detectors and demonstrates that 
measuring the current cross correlations of the two de- 
tectors allows one to get arbitrarily high values of the 
peak-to-background ratio: i.e., it is possible to eliminate 
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the bound that exist in the single-detector case. In Sec- 
tion |Vl we demonstrate how the proposed setup can be 
used to diminish the added position noise of the oscilla- 
tor induced by the presence of the detector, allowing po- 
sition measurement beyond the standard quantum limit 
derived for a single detector. 
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FIG. 1: (Color online) The two typical detector configura- 
tions examined in this article. In both cases, the movement 
of the oscillator is along the x direction in the xy plane, as 
depicted by the ^ sign, a) In-phase configuration, where two 
detectors (with bias Vi and V2, respectively) are located on 
the same side of the central part of the oscillator, such that 
both detectors couple in the same way to the position of the 
oscillator. This is covered in Sec. IIV Al b) Out-of-phase con- 
figuration, where the detectors are located on each side of the 
oscillator. When the position of the oscillator is such that 
the tunneling amplitude of one junction is increased, the tun- 
neling amplitude of the other junction is therefore decreased. 
This is covered in Sec. IIVBI 



II. EQUATION OF MOTION FOR THE 
DENSITY MATRIX 

Approaches based on quantum master equations 
have proven useful in the study of nanomechanical 
systemsSSii^. By writing the equation of motion for the 
density matrix of the full (detector and oscillator) sys- 
tem and tracing out the detector degrees of freedom, one 
can obtain an equation of motion for the reduced density 
matrix describing the evolution of the oscillator taking 
into account the coupling to the detector. In order to 
investigate electronic transport in the coupled system, 
it is useful to refine this approach to keep track of m, 
the number of charges that passed through the detector. 
This allows one to calculate an equation of motion for the 
771-resolved density-matrixJ>2J^i^, a quantum equivalent 
to the m-resolved master equation approach widely used 
in the study of transport properties of classical nanome- 
chanical systems^^i^. 

To study the current cross correlations between two tun- 
nel junction position detectors coupled to an oscillator, 
we use such a fully quantum approach. We label the de- 
tectors with the index a — 1,2 and model each of them 
as a pair of metallic leads with constant density of states 



Aq (in the energy range relevant to tunneling) coupled 
via the tunneling Hamiltonian Htun- The Hamiltonian 
for one detector can therefore be written as a sum of a 
bath Hamiltonian Hs.a describing the leads of junction 
a and a tunneling Hamiltonian i/tun.a 

(1) 

k q 

k,q k,q 

(3) 

where k{q) is a wave-vector in the right(left) lead. The 
coupling between the detector and the position of the 
oscillator is modeled by a linear x— dependence of the 
tunneling amplitude 



Ta{x) 



(to,q 



(4) 



In this equation, to,q. is the bare (oscillator-independent) 
tunneling amplitude of detector a, x is the position op- 
erator of the oscillator and Ti^a is the part of the full 
tunneling amplitude detector a that depends on the po- 
sition of the oscillator. We allow for a general relative 
phase r]a , describing the details of the coupling between 
the tunnel junction and the oscillator. Such a phase can 
in principle be controlled by a magnetic flux penetrating 
an extended tunnel junction consisting of a loop contain- 
ing two junctions, one of which couples to the oscillator. 
Note that in our notation tq^q is dimensionless and ti_q, 
has dimensions of one over length and that we assume 
for simplicity that the tunneling amplitudes do not de- 
pend on the single particle energies Si. ct(q,a)- The opera- 
tor Ya^^ decreases (increases) TOq , the number of charges 
that tunneled through junction a. Its presence in the 
tunneling Hamiltonian allows one to keep track of the 
transport processes that occur during the evolution of 
the system. 

We are interested in calculating the equation of motion 
for the reduced, mc-resolved, density matrix 



p(mi,m2;t) = (toi, r7i2|posc|"ii, 



(5) 



where posc = Trs{ptot} is the reduced density matrix 
that is obtained by tracing out the leads' degrees of free- 
dom from the full system density matrix. Within a Born- 
Markov approximation, the equation of motion of posc 
can be expressed as 

d i 

1 r° 

-J^ dtTlB{[Htun, [Htun{t),Posc{t) <E) Pb]]} , 

(6) 

where Htun = ^^tun,i + ^^tun,2 is the total tunneling 
Hamiltonian, the trace is on both pairs of leads, is 



3 



the coupled density matrix of the two sets of leads and 

= h^ia^a + 1/2) = 1^ + ^^i^ , (7) 
Htunit) = e'^«-*/''i/tu„,oe-^^°-*/'' . (8) 

a 

with Ho^a = ^^osc + Hs^a- In our system, the Born 
approximation corresponds to assuming that tunneling 
in both tunnel junctions is weak enough so that it can 
be treated using second-order perturbation theory. The 
Markov approximation, on the other hand, is valid as 
long as the typical correlation times in the leads (h/eV) 
are much shorter than 27r/r2, i.e. the typical evolution 
time of the oscillator. In practice, this limits the applica- 
bility of the following results to the strongly biased case 
eV 3> hQ. This is experimentally feasible since typical 
oscillator frequencies fl are between 10 - 100 MHz and 
the measurements are done at a much larger bias voltage 
than these frequencies Fii^i^ii^ 

Since the leads of detector 1 are totally independent of 
those of detector 2, can be written as a tensor product 
of the density matrices describing each pair of leads pb = 
PBi <8) Pb2- Also, as Htun.a has no diagonal contribution 
in the basis that diagonalizes Hb^oi, the trace over leads 
q; of a quantity that is linear in i?tun,Q vanishes. As a 
result of those two properties, the trace in Eq. © can be 
rewritten as a sum over two traces, each involving only 
one pair of leads 

TTB{[Htun, [Htunit), Posc{t) (g) Pb]]} 
un,a 1 uii.a 

a 

This effectively makes the two-detector problem two 
single-detector problems. The trace over the leads' de- 
grees of freedom is then carried out in the standard way3£ 
As mentioned above, we are interested in calculating the 
time-evolution of the ma-resolved density matrix. Thus, 
we have to calculate {mi,m2\dtPosc\'mi,m2). We use the 
relations 

{■mi,m2\YiY^ Posc{t)\mi,m2) = p{mi,m2;t) , 
{■mi,m2\Y^ Posc{t)Yi\mi,m2) = p(toi - 1, m2; t) , 
{■mi,m2\Y^Yiposc{t)\mi,m2) = p{mi,m2;t) , 
(ml,TO2|Yl/9osc(^)^l^l"^l,TO2) = p(mi -M,m2;t) , 

as well as the equivalent identities for detector 2 in Eq. ^ 
to find the equation of motion for p(mi, m2; t). 
Following a counting-statistics approac h^^'^'' , it is partic- 
ularly useful to express the equation of motion in terms of 
a counting field Xai the conjugate quantity to the trans- 
fered charge ma- Indeed, Fourier-transforming in the 
transfered-charge indices TOq, 



p{xi,X2;t)= 



m2 = — c 



e'^^"'^p{mi,m2;t) 
(10) 



leads to an equation of motion from which the 
time-dependence of all moments of m (for example, 
dt{'ma),dt{m'^), ...) can be determined. The current- 
current correlations can then be obtained by taking suc- 
cessive derivatives with respect to («%«) of the equation 
of motion of p{xi,X2]t). 

In the regime of weak coupling between the oscillator 
and the detectors, we can write the equation of motion 
of p{xi,X2;t) as 



dt 



P{Xl,X2;t) = [i/osc,p(t)] + ^ ^ [Fai-, p(t)] 



X]7(T,a[i,{p,p(t)}] - [i, [x,p{t)]] 



E 



T-2 
'l.a 



x)p(t){To 



'Tl.aX 



ro,ari,a(e'""°pp(i)-e-'""°p(t)p) 



+ riJppit)^ - xpit)p) 



(11) 



Since p{xi = 0, X2 = 0; t) = Em^ p{mi,m2;t), tak- 
ing Xi = X2 = corresponds to completely tracing out 
the charge degrees of freedom. In this case, one finds 
that /5(0, 0; t) is of Caldeira-Leggett foim^MM., We can 
thus identify the constants Da-^a and jcr.a as, respectively, 
the diffusion and damping constants induced by forward 
((T ~ +) or backward (ct = — ) propagating currents in 
detector a. We can also identify Fa as the average back- 
action force exerted on the oscillator by detector a. We 
find explicitly 



^{—][T,,am)+'^aA-m]. (12) 



Fa 



2Mn \To,a 

i sin(rya 



Tl,Q 



J2 2^^-," 



(13) 
(14) 



where the two inelastic tunneling rates are given by 



(15) 



|2 



(16) 



X I dek,a fk,a{£kA{^^ ^ fq,a{£k,a + E)j , 
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involving a transfer of energy E from the oscillator to 
the lead electron. We denote by T^ ^ the forward tun- 
neling rate, i.e. the rate at which electrons tunnel in the 
direction favored by the voltage bias. The backward rate 
r__Q, corresponds to the reverse process. In Eqs. ([T5|) and 
(fT6|) . fk,a — fR.ai^k.a) is thc Fcrmi distribution function 
describing the local thermal equilibrium of the right lead 
of detector a and fq^a = fh.ai^q.a) is the same for the 
left lead. 

Comparing these relations with the one derived in the 
single-detector casein shows that the full damping and 
diffusion coefficients governing the evolution of the oscil- 
lator are the sum of two single-detector contributions. 
The Caldeira-Leggett form of Eq. pT|) allows us to in- 
clude the effect of direct coupling of the oscillator to the 
environment by adding detector-independent contribu- 
tions Do = 2M7oA:_bTo and 70 = ^/Qo (where Qo is 
the extrinsic quality factor of thc mode) to the previ- 
ously derived diffusion and damping constants. The evo- 
lution of the oscillator is then governed by the two con- 
stants Aot = Do+J2cr,a D-y-.a ^nd 7tot = 7<t,q- 

For the specific case where the electronic temperature 
is zero and where eVa ^ Kl, current will only be pos- 
sible along the {a — -\-) direction, and both 7-, a and 
D^a will be zero. In this case one can also show that 
7+,Q = fi.rf_„/(47rM) and that the diffusion parameters 
are given by Dj^^a = Mj+^asVa- 

The equation of motion for different moments {x^p^) of 
the oscillator can be evaluated by taking the trace of 
p{Q,0]t). More generally, equations of motion for 
combined moments of charge and oscillator quantities can 
be obtained by also considering derivatives with respect 
to the counting fields Xa 

/ «(ii3+n4) 

(17) 



Using the single-detector analogue of Eq. (fTTj) . CG 
showed that, under the conditions mentioned above and 
to first non- vanishing order in ti, the current noise of a 
tunnel junction position detector is given by 



^r(c^) = 2e(/) + ^(2roTi) 



eV 



47r (a;2) 



Sx{uj) 



(18) 



where Aa:;^ = h/{2MQ) is the average of x'^ in the ground 
state of the (quantum) harmonic oscillator and 



Sxiuj) = 



8-ftot^^{x^ 



47 



tot^ 



(19) 



its power spectrum. The full current noise is the sum of 
the usual frequency-independent Poissonian (shot) noise 
and the contribution of interest due to the coupling of 
the junction to thc oscillator. This second part is itself 
expressed as the difference of a classical part (which is 
proportional to y^) and a quantum correction (which is 
proportional to V). 

A relevant figure of merit of such detectors is the peak- 
to-background ratio TZ{uj): the ratio of the contribution 
of the oscillator to the full current noise at frequency 
uj over the unavoidable frequency-independent intrinsic 
detector noise. This ratio is maximal at = and, in 
the case where one only considers the oc contribution 
in Eq. (|18p . was shown to be given by 



5'f*(») - 2e(/) 
MI) 



= 4t.^ 



eV P"^ 
hitot 1 + 



(20) 



where we used (/) = {dtmiit)) ~ e'^Vrg {l + (3'^) /h 
and introduced the dimensionless sensitivity parameter 
= {x'^)/tqM At this point, one should proceed with 
care when maximizing TZ with respect to the sensitiv- 
ity parameter, as 7tot = 7o + 7+ depends on through 
■y^ = {Qt^ /2tt){AxI/ (x'^))l3'^ . Writing out explicitely all 
terms in TZ that depend on /3, one finds that 



III. SINGLE-DETECTOR CASE: BOUND ON 
THE PEAK-TO-BACKGROUND RATIO 

One of the main motivations for studying position mea- 
surements using cross-correlated detector outputs is to 
remove the bound on the peak-to-background ratio that 
appears in the single-detector case, just like in the case 
of a weak measurement of a two- level system^^. In this 
section, we first review the results of Clerk and Girvin 
(CG) fl3| for the single-detector configuration, in the case 
where one considers the dc-biased, T — 0, tunnel junc- 
tion where the a;— dependent tunneling phase is 77 = 0. 
We then carefully analyze the peak-to-background ratio 
and show that this quantity is bounded from above in the 
single-detector case, for finite bias voltage and oscillator 
displacement. 



2tS^ eV 

TT fill 



1 4-/32 



(21) 



is a non-monotonic function of the sensitivity parameter 
(3. For a given (x^), one can then find an optimal value 



27T {x^} 
" Qori Axl ' 



Popt — 



for which TZ is maximal 
_,Qor^feV 



27r V Ml 




Qqt^ Ax l 
2n (x2 



(22) 



(23) 



We can examine this result in two different limits. The 
first is when the damping is mainly detector-independent 
(70 ^ 7-i-)j like in the case where the extrinsic quality 
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factor of the resonator is low, Qo {x^) / {tq IS.Xq) . In 
this case, the maximal peak-to-background ratio, 



measurement^ 



7^ 



eV ( T^Qo Ax, 



Axl Ml 



27r (a;2) J \l + p 



< 4 



Axl Ml 



27r 



is reached when the sensitivity parameter [3 is extremely 
large. However, since the rightmost term of Eq. ([M]) is 
by definition small in this limit, the peak-to-background 
ratio cannot become extremely large when the extrinsic 
resonator damping dominates the detector-induced one. 
Indeed, the real maximum of TZ is reached when one 
considers the opposite limit of a very high resonator Q- 
factor-2, Qo > {^'^) I {'^o ^^o) ■ ^^''^ 7o = the peak-to- 
background ratio can be shown to obey 



eV 1 



Axf^ Ml 1 



< 4- 



eV 



(25) 



In the single-detector case and for given system parame- 
ters {eV and (a;^)), the peak-to-background ratio is there- 
fore always bounded whatever the strength of the cou- 
pling and the bound does not depend on Qo and To. As 
can be seen from Eq. , the peak-to-background ratio 
is in this second case maximal in the limit /3 — > of van- 
ishing coupling. While the optimal TZ can be increased by 
increasing the bias voltage, we stress that our bound on 
TZ denotes the optimal value of the peak-to-background 
reachable for a set of fixed system parameters. 
The nature of the true bound on TZ (i.e., the one found in 
the case Qo ^ cxd) is very similar to the Korotkov-Avcrin 
bound that arises in the context of a weak measurements 
of a qubit. To make this more apparent, we can derive 
this bound following the linear-response approach that 
has been used to derive the bound on TZ in the mea- 
surement of two-levels systems, treating the detector as 
a position-to-current linear amplifier with responsivity 
(dimensionful gain) A ~ 2e^VToTi/h. As noted by CG, 
considering only the dominant cx V'^ term in Eq. (|18p 
corresponds to writing ASi = Sj — 2e(/) = X^Sxii^)- 
At resonance, the power spectrum ASi — 2\^{x'^)/^ 
is inversely proportional to the damping rate 7, in the 
same way that the response of the detector measur- 
ing a qubit is inversely proportional to the dephasing 
rate due to the measurement device. Moreover, in both 
cases one can show that the dephasing (damping rate) 
is proportional to the fluctuations of the bare input of 
the detectors. For a position detector in the high ef- 
fective temperature limit fcsTcff 3> Ml, the detector- 
induced damping is indeed proportional to the sym- 
metrized detector force noisei^ 7 = Sp/^MkBTeS, such 
that ASi < ^MX^{x^)kBTes/'SF- Also, since for a tun- 
nel junction detector there is no reverse gain A' and the 
real part of the cross-correlator Sjp{uj) vanishes-^:, the 
condition on quantum-limited efficiency of the position 





SiSf > — (Rc[A 
4 



X']f + {Re[SiF]y 



(26) 



becomes exactly the one used to derive the Korotkov- 
Averin bound SiSf > We then find that TZ = 

(24) ASi/Si < S{x^)kBT,sl{hnAxl). Using fcsTeff = eV/2 
in the tunnel junction system, this result corresponds 
exactly to Eq. the bound previously derived using 

the equation-of-motion approach. 



IV. PEAK-TO-BACKGROUND RATIO IN 
CURRENT CROSS CORRELATIONS 

Extending ideas from the qubit measurement probleroii, 
we now demonstrate how to eliminate the bound on 
the peak-to-background ratio in a position measurement. 
Calculating the current-current correlations between two 
tunnel-junction position detectors, we show that for cross 
correlation measurements, TZ diverges. We also obtain 
analytical results for the cross correlations in two typical 
cases. 

To calculate the current cross correlations, we use the 
generalized MacDonald formula^i^, a general result 
(valid for stationary processes) that provides a way, in 
the present case, to relate the symmetrized cross correla- 
tions to the Fourier sine-transform of the time-derivative 
of the covariance of mi and m2, the number of charges 
that tunneled through each junction. The generalized 
MacDonald formula reads 



Si,j.,{uj) = 2e^uj dtsin{ujt)Ki^2{t) 



where we defined 
d 
di 



KiAt) 



-{{mim2)t' - {mi)t'{m2)t' 



(27) 



(28) 



In this last equation {'mim2)t corresponds to 
Tr mi m2 p(0, 0, t) and represents the coupled mo- 
ment of TOi and 7712 at time t. 

To proceed further, we restrict ourselves to the case of 
zero electronic temperature and dc-bias. In the following 
subsections, we analyze in detail the two different cases 
depicted schematically in Fig. [1] We have in mind that 
a realization of the setup shown in Fig. [T] is made in a 
similar way as the single-detector setup in Ref. [16. This 
means that the tunnel junctions correspond to atomic 
point contacts (formed by electromigration) which are 
separated by about 1 nm from the oscillator. In contrast, 
the two detectors are assumed to be separated from each 
other by at least 20 nm. Therefore, capacitive cross- 
talking between the detectors will play a negligible rolci^ 



A. In-phase configuration 

We will first consider the case where both 771 = 7^2 = 0, 
the case where both tunnel junctions are located on the 
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same side of the oscillator, cf. Fig. [T]a). To calculate the 
cross correlations, we use Eq. (|17p (with ni — n2 — 0), 
to find that 



dt 



{mim2)t = 
eVi 

-7+,i(m2)t - 7+,2(TOi)t , 

and therefore that Ki 2{t) in this case is given by 

' h 
.eV2 



Ki,2{t) =2^To4Ti,i((a;m2))t + — Ti,i((2;^"i2))t 



+2^ro,2Ti,2((a;mi))f + ^Ti,2((2;^'"i))t • 



closed system of differential equations involving {(pma))t 
and {{xma))t can be generated. This system can be 
solved, using the boundary conditions TOa(O) = and 
assuming that all averages that do not contain TOq are 
time-independent and can therefore be evaluated in the 
stationary [t oo) Iimi1j2^. 

Solving for the different covariances, we find that the 
current cross correlations can be written as 



SfXiio) =e^(2ro,iro^2)(2riari^2) 

'e^ViV2 e{V^+V2) n Axl\ 

7I\ ) '^-(^) ' 



/l2 



A1A2 f 



2h 47r (x2 



nn{eVi+eV2) l^xl 



(31) 
(32) 



(29) 



where the double bracket denotes the covariance of two 
quantities: {{ab))t = {ah)t — {a)t{b)t- This means that, 
to lowest order in ti^q,, the full cross-correlated output of 
the detectors is given in this configuration by 



eVi 



— Ae uj / dtsm{ujt)x 

J72=0 "'0 

eV2 



, TQ^iTi^i{{xm,2))t + —r-To.2Ti,2{{xmi))t 
h h 



(30) 



The cross-correlated signal does not contain any 
oscillator-independent contribution. Using Eq. (|17p . a 



where we introduced the gains Aq, — 
2e^To,aTi,aVaCos{ria)/h. Evidently, the cross-correlated 
output of the detectors (pij) does not contain any 
frequency-independent background noise. The peak-to- 
background ratio Tl{fl) therefore diverges for all values 
of 70/7+7 not because of an increased signal but due to 
the absence of background noise in this configuration. 
For this type of measurement, a relevant figure of merit 
of the detection system TZc is the ratio of the cross- 
correlated output over the frequency-independent noise 
power of individual detectors : TZc — (fi) / V 5*1 6*2, 

where Sa ~ '2e{Ia). For our position detector, we find 



J 



VS1S2 



niTi,2 



-Pi) 
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TI2V, 



e{Vi + V2) (.t2) ^ e{Vi 



V2) {x^ 



Axl 



Ml 



Axl ' 



(33) 



where we used 2xy < (x'^ + y^). From this inequality, 
we see that the maximal cross-correlated output is found 
for (i) twin-detectors (where Ti^i = Ti_2) and (ii) equal 
bias voltages Vi = ¥2- Also, like in the single-detector 
case, TZc is maximal in the limit where there is no ex- 
trinsic oscillator damping 70 and where the correction to 
the average current due to the coupling to the oscillator 
vanishes (/?„ 0). 

Once again it is instructive to compare our value of TZc 
for twin detectors with the equivalent result in the case of 
a weak measurement of a qubit using cross correlationsii. 
In the latter case, the cross-correlated output was shown 
to be limited to 1/2 of the single-detector signal due the 
increased (doubled) detector-induced dephasing. This is 
the same here. 



B. Out-of-phase detection 



We can also analyze the case where one detector couples 
to -I- a; and the other to would happen if the two 

detectors were located on opposite sides of the resonator 
(see Fig. [ij. In terms of the tunneling phases rja, this 
corresponds to taking 771 = and 772 — n. Using Eq. ([?7)l . 
the cross correlations are then given by 



// w ' 

-J^To,iTi^l{{xm2))t J^To^2Tl,2{{xmi))t 



(34) 
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As the coupling between detector 1 and the oscillator is 
the same as in the previous case {{xmi))t remains un- 
changed in this second configuration. The covariance 
{{xm2))t on the other hand changes sign (but keeps the 
same norm) in this new configuration. Equation (j34p 
then yields 



Vl=0 

n2=o 



(35) 



The cross correlations in the second configuration are 
the same as in the first one, but of negative sign. From 
an amplifier point of view, this is easily explained since 
putting r/2 = TT corresponds to transforming A2 —X2 
in Si-^j2 — XiX2Sx- Finally, note that this configura- 
tion was analyzed for two single-electron transistor po- 
sition detectors coupled to a classical oscillator, in Ref. 
(3l| by Rodrigues and Armour. In their article, these 
authors only explicitly calculated zero-frequency cross 
correlations between the currents in both detectors, but 
they conjectured that, at the resonance frequency of the 
oscillator, this detector-configuration (corresponding to 
r]i = 0,r]2 = TT in our approach) should yield strong neg- 
ative cross correlations, just like the ones predicted here. 



BOUND ON THE ADDED DISPLACEMENT 
NOISE 



best possible detector therefore adds exactly as much 
noise as a zero-temperature bath of frequency fii^iM. 
Before discussing the limit on the added displacement 
noise in a cross correlation setup, it is helpful to de- 
scribe how the quantum limit on 5^'^'^(ri) is reached in 
a single-detector configuration. Let's consider for defini- 
tiveness the experimentally relevant configuration where 
eV 3> ksTo > hfl. For a measurement to be quan- 
tum limited, the effective temperature of the oscillator 
Tcs = {j+eV/2 + 7oA;BTb)/(fc_B7tot) must not be dra- 
matically higher than Tq. This is natural, since added 
fluctuations due to the higher effective temperature are, 
by definition, unwanted back-action noise. In this regime, 
one therefore cannot expect S^'^'^ to be close to the quan- 
tum limit unless 7+ <^ 70. The regime of 7+/70 in 
which quantum-limited displacement sensitivity can be 
achieved is therefore very different from the one where the 
bound on the peak-to-background ratio can be reached. 
Using the expression for the full current noise derived 
earlier (Eq. pS]) ). we write the full position noise as 



_ 2e{I) ( hn Ax, 



(37) 



MI) 
A2 



- 8M-/totkBTaa\g{i^)\'^ - 2Af7t, 



' eV 



(38) 



As shown in Sec. IIIIl to derive the equivalent of the 
Korotkov-Averin bound in a position measurement, one 
needs to consider the full current noise, where no dis- 
tinction is made between the signal due to the intrin- 
sic equilibrium ffuctuations of the oscillator S'p-{Ljj) and 
the remainder of the signal S'|^'^(a;). This second contri- 
bution contains, amongst other things, the added signal 
due to heating of the oscillator by the detector. When 
trying to measure precisely the equilibrium ffuctuations 
of a nanomechanical oscillator however, it is important 
to consider the two contributions separately: S^'^{lij) is 
exactly what you would like to measure while S'|'^'^((jj) 
limits the sensitivity of the measurement. When using a 
single linear detector like the tunnel junction, this mea- 
surement sensitivity is quantum-mechanically bounded 
from below^^. 

When discussing this bound on added noise, one usually 
considers the added displacement noise, that corresponds 
to the added current noise referred back to the oscillator. 
We therefore introduce the total displacement noise S] 
defined as 



tot 

X 1 



A2 



^add 



(c) + Sl^{u) 



(36) 



where A is the a;-to-/ gain of the detector, S^'^'^ito) is the 
part of the full displacement spectrum that arises due to 
the presence of the detector. In the relevant limit of a de- 
tector with a high power gain [eV 3> Ml), it was shown 
using general arguments that Sl'^^{^l) > h/Mi}-ftot- the 



where in the last line we introduced the oscillator's re- 
sponse function g~^{u!) = M[{n'^ — w^) -I- 2ijtoti^]- Split- 
ting the second term into a detector dependent and in- 
dependent part, we ffnd 



= 8M7ofcBTo|5(c^)P , 



^add , 



2e{I) 



eV. 



2 +8M7+ — |g(^)|^ 



(39) 

2M7tot^TT^|5(^)l' • 



eV 



(40) 



This way of writing the equilibrium ffuctuations implies 
that we consider 7tot — 7o in 3(^)1 in agreement with our 
previous assumption that 7+ <C 70- The added noise con- 
tains three contributions, corresponding to the detector 
shot noise, the detector-induced heating of the oscilla- 
tor and a correction (oc Hfl/eV) arising from the cross 
correlation between the detector output noise and the 
back-action force, Sjp, respectively. Explicitly, taking 
(/) ~ e'^T^V/h^ we obtain 



add 



eV^ 



hT?eV\g{co)\^^ 



2Af7tot^l5(^)P 



eV 



(41) 



For a ffxed bias voltage, the relevant tunable parameter 
is directly the detector-oscillator coupling ti (and not 
the dimensionless sensitivity parameter (3, since S^'^'^ is 
independent of To)i^ For strong coupling, S^'^'^ is dom- 
inated by heating of the oscillator, while for weak cou- 
pling, the shot noise contribution (cx 1/t^) dominates. 
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This is the regime in which the current generation of ex- 
periments are operated^^. There is an optimal couphng 
''"i.opt = "'/('S^IsC'^)!) that minimizes the total added 
noise. At the resonance, we recover the inequality 



^add 



2eVj 7totM17 



(42) 



where the bound is reached when ri = ri^opt- This is 
the quantum limit on the added displacement noise for 
the single-detector configuration. In passing, we note 
that the effective temperature of the oscillator when the 
coupling strength ti is optimal is 



off 



To- 



m 

4fc^ 



(43) 



in agreement with the general analysis of Ref . [l9| . The 
heating of the oscillator by the detector is, as expected, 
very low when doing a quantum- limited measurement. 
We can now show how cross correlations can be used to 
beat the quantum limit on S^'^'^ derived in the single- 
detector case. In both cross correlation configurations 
(77 — 0, tt), 5*^°* = Si^j^/Xi\2 is identical. Like in the 
single-detector case, we separate the total position fluc- 
tuations in two parts 



5^'i = 8M7ofcsTo|.9(w)P , 



(44) 



add 



= 4 XI 7+.aeK 



7tot - 



{h9.f{eVi+eV2) 



eVieV2 



(45) 



The cross-correlated position spectrum does not contain 
the frequency-independent shot noise contribution that 
diverges for low coupling (c>c l/r^). Therefore, one does 
not need to tune the coupling to equilibrate the "shot 
noise" and back-action "heating" contributions. Instead, 
one can freely tune the couphngs ti^q, such that S'^'^'^{lo) 



vanishes completely. We find S^'^'^ 
47rM7+_Q,opt/^, where 



for 



l,a,opt 



7-|-,a,opt 



7tot 



Ml 



(46) 



At the optimal coupling point, the effective temperature 



of the oscillator is 



Tcff — To 



Ml 

eVi 



Ml 
^2 



Ml 

8fc 



ZT7 + ZT7 (47) 



In the regime of interest {eVa ^ hft), the additional 
heating of the oscillator considerably reduced from the 
single-detector value. 



VI. CONCLUSION 

In this article, we have shown that, for a tunnel-junction 
position detector coupled to a nanomechanical oscillator, 
the optimal peak-to-background ratio TZ at the resonance 
frequency of the oscillator is bounded. In contrast to the 
universal (independent of all system parameters) bound 
derived for a continuous weak measurement of qubits 
{TZ < 4), the new bound derived for position measure- 
ments is a function of the effective temperature of the 
oscillator and its average displacement. We have also 
shown that adding a second detector and using the cross 
correlations between the two detectors allows one to elim- 
inate this bound on TZ. We have analyzed in detail the 
cross-correlated output of the position detectors in two 
typical configurations, and have shown that in both cases 
the optimal cross-correlated signal is measured by twin 
detectors. We also investigated the quantum-limit on the 
added displacement noise and shown that it is possible 
to totally eliminate the added displacement noise by do- 
ing a cross-correlated measurement. This configuration 
therefore opens the door for displacement measurement 
with sensitivities better than the quantum limit. 
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